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Abstract. We consider the break-up of invariant tori in Hamiltonian systems 
with two degrees of freedom with a frequency which belongs to a cubic field. 
We define and construct renormalization-group transformations in order to 
determine the threshold of the break-up of these tori. A first transformation is 
defined from the continued fraction expansion of the frequency, and a second 
one is defined with a fixed frequency vector in a space of Hamiltonians with 
three degrees of freedom. 

1. Introduction. We consider the dynamics of a particle in a potential V de- 
scribed by the following Hamiltonian : 

H{p,x,t)= l -p 2 + V(x,t), (1) 

where V is 27r-periodic in x and t. For V = 0, the phase space is foliated with 
invariant tori. For V sufficiently small and regular, the KAM theorem states that 
most of these invariant tori persist. 

The purpose of this paper is to construct renormalization transformations in 
order to study the persistence of invariant tori with a given cubic frequency to 
smaller than one, satisfying 

u) 3 = auj 2 + bu + l, 

where (a, b) £ 1? . 

For Hamiltonian systems with two and three degrees of freedom, renormalization- 
group transformations in the framework of Ref. |Q have been successfully applied 
for the determination of the threshold of the break-up of invariant tori || ||, ||. 
They have also been used for the analysis of the properties of critical invariant tori 
(at the threshold of the break-up) § § |, % §. 

For two degrees of freedom and for reduced quadratic irrational frequencies, 
renormalization transformations have been defined with a fixed frequency vector, 
i.e. at each iteration of the transformation, the frequency vector of the torus was 
unchanged. This feature was based on the fact that the continued fraction expan- 
sion of the frequency was periodic. For cubic (and more generally for non-quadratic 
irrational) frequencies, this is not the case, and thus ergodic renormalization was 
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constructed from the entries of the continued fraction of the frequency. This trans- 
formation allows the accurate determination of critical couplings by comparison 
with Greene's criterion || and Laskar's frequency map analysis jn| (see Ref. 

In order to define a renormalization with a fixed frequency vector, one has to 
add one degree of freedom to the system as explained below. For degenerate Hamil- 
tonians with three degrees of freedom, a renormalization transformation has been 
defined for the spiral mean frequency vector. This transformation allows one to de- 
termine accurately the thresholds of break-up of tori by comparison with Laskar's 
Frequency Map Analysis || . 

In this article, we propose to define and study numerically two renormalization 
transformations for invariant tori with cubic frequency : one based on the continued 
fraction of the frequency acting within a space of Hamiltonians with two degrees of 
freedom, and the other one with a fixed frequency vector acting within a space of 
degenerate Hamiltonians with three degrees of freedom. 

We apply both renormalization methods to the numerical determination of the 
threshold of invariant tori of a given one-parameter family of Hamiltonians, Es- 
cande's paradigm model Jl^] which is a forced pendulum model. We consider three 
cases : the spiral mean frequency, the Tribonacci mean frequency, and the r-mean 
frequency. 

2. Renormalization from the continued fraction expansion. In order to 
investigate the torus with frequency u>, we map the time-dependent Hamiltonian ([!]) 
into a time-independent Hamiltonian with two degrees of freedom, written in terms 
of actions A = (A\, A 2 ) € M 2 and angles (f — (<fii,f2) £ T 2 : 

H(A 1 ,A 2 ,(p 1 ,(p 2 ) = \^A\ -A 2 + V((px,-(p 2 ). 

The two-dimensional invariant torus is now located around A\ = lo. We shift the 
action Ai in order to locate this torus around Ai — 0. The Hamiltonian becomes : 

H{A, tp) = u ■ A + \^A\ + V{ip x , -ipi), 

where the frequency vector of the considered two dimensional torus is now lo = 
(uj,-1). 

The renormalization we define relies upon the continued fraction expansion of 

LJ = ■ j — = [<z ,ai, . . .]. 

This transformation will act within a space of Hamiltonians H of the form 

H(A,ip) = u-A + V(fl-A,<p), (2) 

where $7 = (l,a) is a vector not parallel to the frequency vector lo. We assume 
that Hamiltonian (|^) satisfies a non-degeneracy condition : If we expand V into 

V(fl-A,<p)= K, (fe) (^' A) k e lv - V 7 

k>0 

we assume that Vq 2 "* is non-zero. 

The transformation contains essentially two parts : a rescaling and an elimina- 
tion procedure. We follow the renormalization scheme defined in Ref. [pi. Similar 
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renormalization transformations have been defined in Rcfs. [[l3| |TJ| for vector fields. 

(1) Rescaling : The first part of the transformation is composed by a shift of 
the resonances, a rescaling of time and a rescaling of the actions. It acts on a 
Hamiltonian H as H' — H o T (see Ref. ^ for details) : 

H'(A,<p) = Xoj-'H f-iiV A,-JV-Vj , (3) 

with 

A = 2co- 1 (a + a) 2 V (2) , (4) 

and 

N a =(l J), (5) 

and a is the integer part of uj^ 1 (the first entry in the continued fraction expansion). 
This change of coordinates is a generalized (far from identity) canonical transfor- 
mation and the rescaling A is chosen in order to ensure a normalization condition 
(the quadratic term in the actions has a mean value equal to 1/2). For H given by 
Eq. (||) , this expression becomes 

H'(A, tp) = J ■ A + K [k) («' ' A) V 1 "^ (6) 



v.k 



where 



u>' = (w',-1) with J = oj- 1 - a, (7) 
SI' = (l,a) with a' = l/(a + a), (8) 

= r fc ^l with r fe = (-l) fc 2^ V' 2 (a + a) 2 ~ fc (Vo (2) ) . (9) 

We notice that the frequency of the torus is changed according to the Gauss 
map : 

LJ H-> J = UJ- 1 - [iJ- 1 ] , (10) 

where denotes the integer part of a; -1 . Expressed in terms of the continued 

fraction expansion, it corresponds to a shift to the left of the entries 

lo = [a ,ai,a 2l ■ ■ ■} i-> u>' = [oi, a 2 , a 3 , . . . ]. 

Remark: After renormalization (one or more iterations) the frequency vector 
w = (u>, — 1) is changed into u/ = fj,Mu) with M G GLi^E). In order to define a 
renormalization with a fixed frequency vector uj' = uj, the frequency w has to be 
a quadratic irrational. For a cubic frequency (e.g., satisfying uj 3 — au; 2 + bui + 1), 
such renormalization cannot be defined with a fixed frequency vector. The idea is 
to add one more dimension and consider frequency vectors of the form (u), — 1, Wi). 
If cl>i is carefully chosen, there exists /i and M € GL^^L) such that a; = fiMu, and 
then one can define a similar renormalization with the exception that now there is 
an additional degree of freedom corresponding to the additional frequency u\. 

(2) Elimination : The second step is a (connected to identity) canonical trans- 
formation U that eliminates the non-resonant modes of the perturbation in H' . 
Following Ref. ||, we consider the set I~ C Z 2 of non-resonant modes as the set 
of integer vectors v = (^1,^2) such that \i/ 2 \ > Other choices of resonant 

modes can be done without affecting the convergence and speed of the algorithm 
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(see Ref. [[TI] for a detailed discussion). A mode which is not an element of I~ , will 
be called resonant. The canonical transformation U is such that H" = H' oU does 
not have any non-resonant mode, i.e. it is defined by the following equation : 

r(H'oU)=0, (11) 

where I - is the projection operator on the non-resonant modes; it acts on a Hamil- 
tonian (||) as : 

l~H= V„ k) {n- Afe lv -*. 

fe>o 

We solve Eq. ([ll]) by a Newton's method following Ref. g [| . We notice that in 
Refs. j^, H a renormalization procedure acting within a space of quadratic Hamil- 
tonians in the actions has been constructed. This renormalization is the one used 
in this paper for numerical purposes. This transformation gave the same accuracy 
than the renormalization defined for Hamiltonians in power series in the actions 
for the numerical computation of the parameters that characterize the break-up of 
invariant tori. 



These rescaling and elimination procedures define the renormalization trans- 
formation acting on a space of Hamiltonians (||) with two degrees of freedom as 
H" = TZ(H) = H o T oU. The main point is that the frequency u) is changed 
according to the Gauss map (|l^) each time we iterate the transformation. 

3. Renormalization with a fixed frequency vector. In order to study a torus 
with frequency lo, we consider another frequency lo\ such that (l,u>, u>i) is a com- 
plete integral basis of the cubic field to which uj belongs to. We consider the invari- 
ant torus with frequency vector uj — (uj,— 1,u>i). The Hamiltonian (Q) is mapped 
into the time independent Hamiltonian with three degrees of freedom written in 
terms of actions A — (A x , A 2 , A 3 ) e E 3 and angles <p — (^1,^2,^3) £ T 3 : 

H{A 1 ,A 2 , A 3 , tpx, tp 2) v? 3 ) = \a\ -A 2 + V((p u -(f 2 ). 

We shift the action A\ by u> and we add W1A3 since A${t) is constant (H does 
not depend on 993). The Hamiltonian becomes : 

H(A, ip) = uj ■ A + -A\ + V(ipi, -ip 2 ). 

We assume that uj is an eigenvector of a matrix N with integer coefficients, 
determinant ±1, and a real eigenvalue /1 of modulus smaller than one. For instance, 

if LOi = bJ 2 , 

/ 1 

N= -1 

\ b -1 a 

and /i = lo. 

This transformation will act within a space of Hamiltonians H of the form 

H(A,<p) = w- A + V(il- A,<p), (12) 

where f2 € M 3 is a vector of Euclidean norm one (||n|| = 1), not parallel to the 

frequency vector us. We assume that Hamiltonian (O) satisfies a non-degeneracy 
(2) 

condition (Vq is non-zero). 

The renormalization is defined by the rescaling and elimination steps. 
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(1) Rescaling : The first part of the transformation is composed by a shift of 
the resonances, a rescaling of time and a rescaling of the actions. It acts on a 
Hamiltonian H as H 1 = H o T : 

H'(A, tp) = X^H [\NA, N- Vj , (13) 

where N denotes the transposed matrix of N, and A is the rescaling coefficient : 

A = 2 M - 1 ||iVO|| 2 V/ (2) . (14) 
For H given by Eq. (tL2h, this expression becomes 

H'(A, tp) = u ■ A + K {k) (n' • A) k e lv - v , (15) 



'.k 



where 



Vi^ = r h V$l with r k = 2 1 -V- 2 ||A^|| 2 - fe (vf 5 )* * . (17) 

We notice that the frequency vector uj of the torus is kept fixed by the transfor- 
mation since from Eq. (jl^), we have A/i _1 u; • \NA = uj ■ A. 

(2) Elimination : The second step is a canonical transformation U that elimi- 
nates the non-resonant modes of the perturbation in H' . 

Following Ref . (J| , we consider the set I~ of non-resonant modes to be the set of 
integer vectors v £ Z 3 such that \u>-u\ > ^||w||||i/||. The canonical transformation 
U is such that H" = H' oU docs not have any non-resonant mode, i.e. it is defined 
by the following equation : 

l-{H'oU)=Q. (18) 

We solve Eq. ( |l8| ) by a Newton's method following Ref. |3j, and we apply the 
procedure acting on a space of Hamiltonians quadratic in the actions. 

These two procedures define the renormalization transformation acting on a 
space of Hamiltonians ( |l2| ) with three degrees of freedom as H" = 1ZH = HoToU. 
The main point is that in this renormalization the frequency vector uj is kept fixed 
at each iteration of the transformation. 

4. Numerical results. We consider Escande's paradigm model jl^] which is a 
forced pendulum model : 

H e (p, x, t) = —p 2 — e(cosx + cos(x — tj). (19) 

We map this Hamiltonian into a time-independent Hamiltonian with two degrees 
of freedom of the form (^|), or into a time- independent Hamiltonian with three 
degrees of freedom of the form (|l2|). We compute the critical coupling at which 
the torus is broken by the two renormalization methods 1Z. We determine e c by 
iterating the renormalization transformation on H e . If the coupling e is smaller 
than a critical value e c , the iterations converge to some integrable Hamiltonian Hq 
(the coefficients V v for v ^ vanish), and if e is larger than e c , the iterations 
diverge : 

K n H £ -> H (A) for |e| < e c , 
Tl n H F -> oo for lei > e,, 
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as n tends to infinity. For numerical convenience, each renormalization 1Z is de- 
fined for quadratic Hamiltonians (all the iterations TZ n H E are quadratic in the 
actions), and at each step, we truncate the Fourier series by neglecting all the 
modes v = (fj)j=i,2 or 1=1,2,3 such that maxi \vi\ > L. Other way of truncating 
the Fourier series can be chosen. The one chosen here is the simplest choice for the 
numerical implementation. 

Since the renormalization for three degrees of freedom involves more Fourier 
coefficients (proportional to L 3 , compared to L 2 for the renormalization for two de- 
grees of freedom), we expect the resulting transformation to be much slower than 
the renormalization for two degrees of freedom. However, it turns out from numer- 
ical computations that in order to reach a given accuracy, both methods require 
approximately the same amount of computational time, i.e. the cut-off parameter 
L for three degrees of freedom is chosen smaller than the one for two degrees of 
freedom. In what follows, the thresholds are computed with approximately the 
same amount of computational time. 

We apply the above procedures for three different frequencies : the spiral mean, 
the Tribonacci mean, and the r-mean frequencies. 

4.1. spiral mean torus. We consider an invariant torus with frequency u> = a^ 1 
where a is the spiral mean i.e. it is the real solution of a 3 = a + 1. It is approx- 
imately a s» 1.32472. The frequency cr _1 belongs to the complex cubic field with 
negative discriminant —23. The break-up of KAM tori with spiral mean frequency 
of Hamiltonian systems with two degrees of freedom has been previously studied 
using scaling law analysis of Greene's residues of nearby elliptic periodic orbits |]l6| . 

Renormalization for u) = (ct -1 ,— 1). The first entries of the continued fraction 
expansion of cr" 1 are [1,3, 12, 1, 1,3,2,3,2,4,2, 141,80,2,5, . . .]. For a cut-off pa- 
rameter L = 25, the renormalization defined in Sec. || gives s c w 0.0155. 

Renormalization for u> = (<x _1 , — 1, c~ 2 ). The frequency vector u> = (cr -1 ,— l,c ) 
is an eigenvector of the matrix 



For a cut-off parameter L = 15, applying the renormalization defined in Sec. |^ 
with the above matrix N, we obtain e c « 0.016. If we choose u)\ = a~ 2 + o~ l + 1 
instead of u>i = o~ 2 , the renormalization defined by the matrix 



gives the same result (e c « 0.016) as the renormalization defined for u>\ = a~ 2 . 

Up to the numerical precision, the renormalizations for two and three degrees of 
freedom give the same result for the critical threshold e c . 

4.2. Tribonacci torus. We consider an invariant torus with frequency 10 = <r _1 
where s is called the Tribonacci number and is the real solution of c 3 = ^ 2 + <r + 1 . 
Its value is approximately ? « 1.83928. The frequency <; _1 belongs to the complex 
cubic field with discriminant —31. 
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Renormalization for u: = (? , — 1). The first entries of the continued fraction ex- 
pansion of s' 1 are [1, 1, 5, 4, 2, 305, 1, 8, 2, 1, 4, 6, 17, 5, 1, . . .]. For a cut-off parame- 
ter L = 25, we obtain e c « 0.02186. 

Renormalization for ui = ,—l,<; ), The frequency vector w = (? , — 1,? ) 
is an eigenvector of the matrix 



For a cut-off parameter L = 10, the renormalization with the above matrix N 
gives e c w 0.0219. We notice that in this case, a very good agreement between both 
renormalizations is obtained even for small cut-off parameters L. 

4.3. r-mean torus. We consider an invariant torus with frequency u> — r _1 where 
t is the real solution of r 3 = -r 2 + 2r + 1. Its value is r = 2 cos(27r/7) w 1.24698. 
The frequency r _1 belongs to the real cubic field with positive discriminant 49. 

Renormalization for u) = (t -1 , —1). The first entries of the continued fraction ex- 
pansion of t -1 are [1, 4, 20, 2, 3, 1, 6, 10, 5, 2, 2, 1, 2, 2, 1, . . .]. For a cut-off parameter 
L = 25, we obtain e c w 0.01247. 

Renormalization for oj = (t , — 1,t ). The frequency vector u; = (t , — 1,t ) 
is an eigenvector of the matrix 



with eigenvalue /i = t _1 . For any cut-off parameter, the critical value we obtain 
by the renormalization procedure for this matrix N (e c ~ 0.48) is very far from the 
one obtained by the previous renormalization. The reason is that the matrix N 
has not the good properties. In order to define the procedure, the renormalization 
transformation has to contract all the integer vectors in / + (the complement of 
I~ in Z 3 ) into vectors in I~ [Q. The spectrum of N is composed by r _1 , —(1 + 
t -1 ) -1 , — (r + 1). The shift of the resonances is a map acting on integer vectors as 
v i — ► A r_1 i/. Since there are two eigenvalues of modulus larger than one for TV -1 (r 
and — 1 — r _1 ), there are two unstable directions for the shift of the resonances. 

We notice that for the negative discriminant case, this problem cannot occur. If 
there is one eigenvalue of modulus smaller than one, then the two other complex 
conjugated eigenvalues are of modulus greater than one since the determinant of N 
is ±1. 

Renormalization for ui — (r _1 , — 1, r + 1). The frequency vector u) = (r _1 , — 1, r + 
1) is an eigenvector of the matrix 



with eigenvalue \i = (1 + t)" 1 . Here the matrix N has the good property of having 
only one eigenvalue of modulus smaller than one (the two other eigenvalues are — r 
and 1 + r _1 ). The renormalization defined in Sec. |3[ for a cut-off parameter L = 7, 
gives e c w 0.013. For this frequency vector, an approximate renormalization has 
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been set up for degenerate Hamiltonians with three degrees of freedom in Ref. [|j 
in order to investigate self-similar properties of three-dimensional critical tori. It is 
defined from the above matrix N and structurally stable dynamics has been found 
for this approximate renormalization. 

Remark : For the forced pendulum model ( |l9| ) , the Tribonacci torus is the most 
robust between the three invariant tori investigated in this paper, and the spiral 
mean torus is more robust than the T-mean torus. However, this feature depends on 
the perturbation in a way that is not understood. The most robust invariant tori are 
conjectured to be noble tori in general (see Ref. Q for area-preserving maps and 
Rcf. for the stochastic ionization of Hydrogen driven by microwaves). Therefore, 
tori with cubic frequency should be less important for the large-scale diffusion of 
trajectories. However, the results presented in this article would be useful in order 
to investigate the break-up of invariant tori in Hamiltonian with three degrees of 
freedom and its comparison with the break-up in two degrees of freedom. 

5. Conclusion. We have constructed and studied numerically two renormalization 
transformations in order to compute thresholds of break-up of invariant tori with 
a cubic frequency in Hamiltonian systems with two degrees of freedom. The main 
advantage of the procedure acting in a space of degenerate Hamiltonians with three 
degrees of freedom, is that it is independent of the continued fraction algorithm 
and it is defined for a fixed frequency vector. Both methods allow an accurate 
computation of critical couplings. 
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